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A vicious walker system consists of N random walkers on a line with any two walkers annihilating 
each other upon meeting. We study a system of N vicious accelerating walkers with the velocity 
undergoing Gaussian fluctuations, as opposed to the position. We numerically compute the survival 
probability exponent, a, for this system, which characterizes the probability for any two walkers 
not to meet. For example, for N — 3, a — 0.71 ± 0.01. Based on our numerical data, we conjecture 
that ^N(N — 1) is an upper bound on a. We also numerically study N vicious Levy flights and 
find, for instance, for N = 3 and a Levy index fi = 1 that a — 1.31 ± 0.03. Vicious accelerating 
walkers relate to no-crossing configurations of semiflexible polymer brushes and may prove relevant 
for a non-Markovian extension of Dyson's Brownian motion model. 

I. INTRODUCTION 

Consider N random walkers on a line. The only interaction between the walkers is that if any two walkers meet, they 
annihilate each other, hence these walkers are deemed vicious. This system was introduced by M. E. Fisher in his 1983 
Boltzmann Medal lecture and was applied to interfacial wetting in 1 + 1 dimensions since the interaction of different 
interfaces (walkers) drives the wetting transition 1, 2] . In addition, the fermionic nature of vicious walkers provides 
a Coulomb gas description and, thereby, a link with Gaussian random matrices 0, 0|. More specifically, Baik has 
proven that a particular limiting conditional distribution of the displacement of the leftmost walker is equivalent to the 
Tracy- Widom distribution for the Gaussian Orthogonal Ensemble(GOE) [5j. Moreover, vicious walker configurations 
correspond to directed polymer brushes with the vicious mechanism capturing the non-intersecting property of the 
polymers [6[ . So while the study of vicious walkers has attracted attention from the mathematical physics community, 
a number of different physical applications also drive its study. 

There have been generalizations of vicious walkers to dissimilar walkers [7], walkers with drift [8], and walkers 
with external potentials [9]. Here, we introduce a system of N accelerating vicious walkers and study the survival 
probability, s(t), the probability of having none of the N walkers annihilated up to time t. In one-dimension, a 
randomly accelerating walker, x(t), is defined by 

d 2 x(t) 

where r](t) is Gaussian noise with (f](t)) — and (rj(t)rj(t')) — 2DS(t — t'). The corresponding Fokker-Planck equation 
is 
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p(x,v,t) = — p(x,v,t), (2) 



where p(x,v,t) is the probability density distribution for a randomly accelerating walker in one-dimension. 

A randomly accelerating walker is presumably the simplest non-Markovian stochastic process. Given the fermionic 
interaction between vicious accelerating walkers, one can explore the possibility of a non-Markovian analogue to 
the Coulomb gas and, hence, a potentially new class of non-Markovian random matrices. Recently, Fukushima 
and collaborators have revisited Dyson's original Brownian motion model for random matrices and found another 
non-Markovian stochastic process after generalizing the coefficients of the matrices [ioj . In addition, a randomly 
accelerating walker appears in the Boltzmann weight of an extensible, semiflexible polymer of length L with non-zero 
bending energy [TT1 |. For a given displacement vector r(z) from some reference point along a contour length, z, the 
Hamiltonian, H , is given by 

H(r(z), u(z); z) = J t dz (^f] , (3) 



2 Jo V dz 

where u(z) denotes the tangent vector and n characterizes the bending rigidity. Mapping the contour length z to time 
t and the displacement vector r to x(t), the equation of motion for this system corresponds to a randomly accelerating 
walker in the corresponding dimension. Implementing the vicious interaction between N random accelerating walkers, 
therefore, corresponds to the statistics of non-intersecting semiflexible polymer brushes. 
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II. N < 3 VICIOUS ACCELERATING WALKERS 



We begin by considering N — 2 vicious accelerating walkers in one-dimension with equal diffusion constants D\ = 
D 2 = D. The two accelerating walkers are governed by the equation, 



d 2 d 2 \ ( d d \ d 

D ( q^2 + q^2 J p(xi,x 2 ,v 1 ,v 2 ,t) - y v i-Q^- +V2 ~Q^~J P( x i, x 2,v 1 ,v 2 ,t) = —p(x 1 ,x 2 ,v 1 ,v 2 ,t) (4) 



with the initial condition, p(x±,x 2 ,vi,v 2 ,t = 0) = S(xi — x\ t i)8(x 2 — x 2t i)S(vi)S(v 2 ), with x\j < x 2 ^. In addition, to 
compute the survival probability, s(t), we implement the boundary condition, 

p(xi = x 2 , vi > i>2, t) = (5) 

such that the system "dies" when the two walkers meet in space and have ingoing velocities. Note that the boundary 
condition in v is redundant since, due to the initial positions, the two walkers cannot meet with a relative outgoing 
velocity. 

We then choose the relative coordinate system to reduce the N — 2 vicious accelerating walkers to one random 
accelerating walker in the presence of an absorbing wall. This change of variables takes the form, x = x\ — x 2 and 
v — v\ — v 2 . The Fokker-Planck equation thereby reduces to 
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dv 2 dx 



d 

p[x,v,t) = — p(x,v,t), (6) 



with the boundary condition, p(x = 0, v > 0, t) = 0. 

The survival probability for this process is nontrivial in that one cannot invoke the method of images as is done 
for the ordinary random walker. Using properties of the integral of a Brownian curve, Sinai (l2l | proved that the 
asymptotic survival probability distribution is given by 

a(t) ~ (7) 

at long times. Note that the first-passage time distribution, f(t), where the first passage time is defined by the time 
at which any of the two walkers meet is given by f(t) = —^p- such that f(t) ~ i~ 5 / 4 . In general, if the survival 
probability distribution is given by s(t) ~ t~ a at large times, then f(t) ~ with /J = a + 1. 

A heuristic argument, based on Sinai's approach, for a = 1/4 was given in Ref. First, a new time counter, 

M, is defined by each "original" time the velocity crosses zero. In other words, the original time is now defined by 
the distribution of first-passage time of a random walk undergoing a Levy flight with Levy index /i = 1/2. Moreover, 
with this new counting, the position variable is also a Levy flight with Levy index fi = 1/3. By invoking the powerful 
superuniversality of the Sparre- Andersen theorem in one-dimension 6j j , the first passage time distribution for the 
position in terms of M is the same as for a random walk, i.e. M~ 3 / 2 . To convert back to the original time, one simply 
needs to compute the integral 

/(*)~ / -^^eM~t 2 /M)dM, (8) 
J M 2 t 2 

where -^f- exp(— t 2 /M) is the limiting distribution for the sum of M fi = 1/2 Levy variables. We should also mention 

that Burkhardt [ll| made the use of Marshall- Watson functions to solve for the Laplace transform version of Eq. 6 
with the absorbing boundary condition. 

To numerically check for a = 1/4 result (and other results) , we implement the go-with-the winners algorithm [T3] ■ 
We do this because as N increases, the first passage time exponent j3 increases, making it more difficult to sample the 
tail of the distribution. The go-with-the winners algorithm iterates replica systems in parallel. Once at least one pair 
of random accelerators have crossed paths in some fraction of the replicas, those surviving replicas are copied over to 
the replicas that have already "died" . We choose that fraction to be one-half such that each copy generated carries 
a relative weight of 1/2 C , where c is the number of copies, which is then incorporated into the survival probability. 
The number of replicas range between 1,000 and 10,000. The number of runs averaged over range between 10 and 
40. Fluctuations between the runs are then used for error analysis. 

Having calibrated our results for the N = 2 case, we now consider N = 3 vicious accelerating walkers with equal 
diffusion constants. The Fokker-Planck equation is given by 
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p(X,V,t) = -^p(X,V,t), (9) 
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FIG. 1: Log- log plot of survival probability distribution versus time for TV = 2 and N = 3 vicious accelerating walkers in line 
and for a single accelerating walker in a 60° wedge geometry. The line denotes a survival probability exponent of 1/4, while 
the dashed line denotes a survival probability exponent of 3/4. 



where X = {xi,x 2 ,x 3 }, and V = {vi,v 2 ,v 3 }. 

We apply change of variables v = v± — v 2 , u — v 2 — 1)3, x = x\ — x 2 , y — x 2 — x 3 such that the LHS operator for 
two relative coordinates becomes 
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\dv 2 du 2 dvdu J \ dx dy J 



(10) 



with the absorbing boundary at x = 0, y = 0, u > and v > 0. Again, here, boundary conditions on the velocities 
are redundant. To remove the coupled term in u and v, we perform another set of linear transforms, 

u = t^l^ = -l-q/V3 ^ x = -z-w/V3 ^ and y = z-w[V3 tQ obtain 

/ d 2 d 2 \ ( d , d \ 

6D {W + 8P)-{ q 8^ + l d- z ) (U) 

with absorbing boundaries, z — ±iy/v3, i.e. a 60° wedge in the z — w plane. 

We have reduced three vicious accelerating walkers in one-dimension to one accelerating walker in two-dimensions 
in a wedge geometry. For comparative purposes, let us review the ordinary random walker in a wedge [i~8j . Using 
the conformal mapping z 1 = z 1 ' 9 , one can map the wedge geometry to the upper-half plane. Then, the motion in 
x-direction becomes unbounded and in the y-direction it is the simple situation of a one-dimensional walker with 
an absorbing boundary condition. The survival probability distribution asymptotically in the upper-half plane is 
s(t) ~ t~ x l 2 . After an inverse conformal mapping, one arrives at the survival probability distribution of s(t) ~ t^^l 26 . 

While there is currently no analytical solution for the survival probability distribution for an accelerating walker in 
a 60° wedge geometry, for a 90° wedge, the Sparre- Andersen theorem can be invoked for the two indepedent directions 
to arrive at an asymptotic power-law survival probability distribution with a = h. In addition, for a 180° wedge, 
a = j. While the Sparre- Andersen theorem is quite powerful and can easily be extended to as many independent 
dimensions as needed, the 60° wedge geometry couples the two directions (for N = 3, at least) and, therefore, the 
Sparre- Andersen theorem, as it stands, cannot be invoked. 

However, to smoothly interpolate between the 90° and 180° cases, we conjecture that for other wedge angles, 
the survival probability distribution also asymptotes to a power-law with survival probability exponent, a, with a 
decreasing continuously as the wedge angle increases. To this conjecture, we resort to numerical simulation of both 
the wedge geometry for one accelerating walker and the line geometry for three accelerating walkers. The result is 
presented in Fig. 1. We measure a survival probability exponent of a = 0.71 ± 0.01 for N = 3 vicious accelerating 
walkers in a line, which agrees with the 60° wedge geometry (with essentially the same error bar). Also, referring back 
to the accelerating walker-Levy flight mapping implemented to demonstrate the a = j exponent for an absorbing 
accelerating walker in one-dimension, simulating a /1 = 1/3 Levy flight in a 60° wedge geometry with the opening 
angle, 6, between 0° and 60°, yields a = 0.71 ± 0.03. 

Figure 2 tests our conjecture that a continuously decreases with increasing wedge angle. Both the numerical values 
of a for the 90° and 180° wedges agree well with their analytical counterparts. For comparison purposes, we have also 
plotted the curve a — ir/46, which agrees with the two analytical solutions and can be viewed as a trivial extension 
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FIG. 2: Survival probability exponent for a two dimensional accelerating walker in a wedge geometry as a function of the 
opening angle. The line denotes a — Inset: Log-log plot to emphasize the difference between the data and a = The 
error bars are smaller than the symbols. 



of the random walker solution. While the agreement looks reasonable at larger angles, the deviation becomes more 
apparent at smaller angles. It also appears that the divergence in a as 9 decreases to zero is slower than 1/6. 



III. N < 3 MIXED CASE 



Before moving on to N > 3, there exists another problem that may be of interest to those studying mixtures of 
flexible and semiflcxible polymer brushes and where some analytical results can be obtained. Consider a mixture of 
vicious accelerating and random walkers — an example of a vicious walkers system with dissimilar members. For one 
random walker and one random accelerator, the random walker can be approximated as an absorbing boundary at rest 
since the displacement of a random walker scales as x(t) ~ t 1 / 2 , while for a random accelerating walker, the velocity 
scales as v(t) ~ t 1 ^ 2 such that the displacement scales as x(t) ~ t 3 ^ 2 . Thus, in long time limit, the displacement of 
a normal random walker is negligible compared with that of a random accelerator. In other words, one can simply 
take the normal random walker as a fixed absorbing wall to fulfill the vicious mechanism. Simulations of one random 
walker and one accelerating walker that annihilate upon crossing verify this. 

Now to address the N — 3 mixed case. We use the "A" to denote an accelerating walker and "R" to denote a normal 
one in the present discussion. The two pure states "RRR" and "AAA" are trivial with the former following Fisher's 
general solution and latter demonstrated by our numerical simulation described earlier in the paper. The combination 
with two normal walkers and one accelerator can have two forms, namely "RRA" and "RAR". The former can be 
treated as one absorbing random walker and one absorbing accelerating acting independently such that a = 3/4. 
The latter is an accelerating walker sandwiched between two random walkers, i.e. a bounded domain, and hence the 
survival probability decays instead exponentially in time. A detailed analysis can be found in Ref.[19(. Simulations 
support both results. See Figure 3. For "ARA", there is an absorbing wall in between two accelerating walkers, i.e. 
a pair of decoupled accelerator-wall systems, hence a = 1/2. The last possible configuration, "AAR", is nontrivial in 
that it can be reduced to one accelerating walker in a 45° wedge geometry, just as two vicious random walkers in the 
presence of an absorbing wall can be mapped to one random walker in a 45° wedge geometry. We measure a nontrivial 
survival probability exponent of a = 0.89 ± 0.01 for this case, which is consistent with our wedge measurements. 



IV. TV > 3 VICIOUS ACCELERATING WALKERS 



Now, we numerically address N > 3 vicious accelerating walkers. To compare with ordinary vicious walkers, 
based on the method of images, Fisher [l], 0| considered one compound walker in N dimensions that cannot cross 
any of the x\ = X2,X2 = Xs,..., or xjv-i = xn linear manifolds. Using the method of images, as long as the 
initial positive and negative weights (corresponding to unrestricted positive and negative walkers) are chosen such 
that probability distribution is antisymmetric under reflection within each linear manifold, then the distribution will 
satisfy the absorbing boundary conditions for all future times. These weights can be represented as a Vandemond 
determinant, which can be factorized to yield a product of ^N(N — 1) pairings such that a(N) = N(N — l)/4, where 
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FIG. 3: Log-log plot of the survival probability distribution for the N = 3 mixed cases of vicious accelerating walkers (A) and 
vicious random walkers (R). The line denotes a first passage time probability exponent of f3 = |. 

a(2) = 1/2. Previous and the current simulations verify this prediction. See Figure 4. By the same argument, if 
the survival exponent for two (one pair of) vicious accelerating walkers is 1/4, the survival exponent for N vicious 
accelerating walkers system should be N(N — l)/8. However, for vicious accelerating walkers, the method of images 
fails. 

Based on the N = 2 and N — 3 results, combined with the fact that N vicious accelerating walkers can be mapped 
to one accelerating walker in N — 1 dimensions in an unbounded domain, we expect the power-law survival probability 
distribution extends to N > 3. In Figure 4 we present the simulation results of survival probability exponents for 
vicious accelerating walkers system up to N — 10, as well as the too trivial prediction, a — N(N~ l)/8. The deviation 
is apparent at large N. We find that N(N — l)/8 is clearly an upper bound for the measured N. The non-Markovian 
nature of the accelerating walker enables the system to survive longer and in a way that cannot be accounted for by 
individual pairings. 
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FIG. 4: Survival probability exponent a for vicious accelerating walkers (solid squares) and vicious Gaussian walkers (solid 
circles) systems up to iV = 10. The simulation results for vicious Gaussian walkers agree very well with theory prediction 
a = N(N — 1)/4 (black curve). However, the results for vicious accelerating walkers deviate from a method of images prediction 
of a = N(N - l)/8 (blue dashed curve). 
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FIG. 5: Log-log plot of the the survival probability distribution for N = 3 vicious Levy flights for different values of the Levy 
index, fi. The curve denotes a survival probability exponent of a = |. Inset: The survival probability exponent for N = 2 and 
N = 3 vicious Levy flights as a function of the Levy index. For /i > 2, we obtain the vicious random walker results. 



V. VICIOUS LEVY FLIGHTS 



As mentioned previously, vicious Gaussian walkers problem is closely related to the Gaussian random matrix 
theory. A matrix with random entries falls into this category as long as the entries are independent and identically 
distributed (iid) variables with a finite second moment of the corresponding distribution. A generalization of the 
random Gaussian matrix is the random Levy matrix .20], where the entries are drawn from a broader distribution, 
namely a Levy distribution. The most important characteristic of Levy distribution is a heavy power-law tail step- 
size S distribution, P(S) ~ S I_1_M for large S. When ji > 2, the variance of the distribution is finite and, hence, 
the central limit theorem holds for the distribution of the sum of independently drawn Levy variables. Similarly, 
the random Levy matrices reduce to random Gaussian matrices. However in the regime (i < 2, the variance of the 
distribution diverges and hence random Levy matrices behave qualitatively different from random Gaussian matrices. 
For example, the famous Wigner-Dyson semicircular law is replaced with a density of states that extends over the 
entire eigenvalue axis (2fj |. 

Inspired by the connection between vicious Gaussian walkers and random Gaussian matrices as well as the connection 
between Levy flights and random accelerating walkers, we study the problem of vicious Levy flights. Given N Levy 
flights in one-dimension, we define the vicious interaction between pairs. Because Levy flights are nonlocal, two Levy 
flights jump over each other without meeting at some exact point. Hence, there could be two ways to define the vicious 
interaction, to prohibit the jump-overs or to allow for jump-overs and the annihilation occurs upon intersection within 



some range, irrespective of the ordering. The latter has been recently studied[2l|. However, we are more interested 
in the former case, in which the set of Levy flights annihilate whenever a crossing occurs. In other words, a surviving 
system strictly maintains the initial ordering of all flights, which is the same as for the vicious Gaussian walkers. 
The Fokker-Planck equation for a system of N one-dimensional vicious Levy flights is described as 

N f)^ Ft 

Ij8SF p(Jt ' 1) = s«"' < 12) 



where the normal Laplacian is replaced by the Riesz-Feller derivative of fractional order 2 > \i > [22J, |23| . This 
derivative has an integral representation, which more easily reveals its nonlocal nature. The initial condition is still 
p(X,t = 0) = rij=i 3{ x ~~ x j,i)> w ith < %k,i f° r au 3 < The boundary condition for the non-crossing vicious 
interaction as we described above is then p(%j, t) = 0, if Xj > Xk for any j < k. 

The N = 2 case is, again, equivalent to the first-passage problem of a single Levy flight via a transformation to 
relative coordinates (and integrating out the center of mass coordinate). The only difference with a random walker 
is that the absorbing boundary condition at the origin has to be modified to an absorbing region occuping the 
positive x-axis to preserve the non-crossing property. The first-passage property of a Levy flight is governed by the 
Sparre- Andersen theorem (l34l6j ■ which implies that the first- passage time distribution for any symmetric step size 
distribution in one-dimension asymptotes to the same as that of a Gaussian walker. Thus, the survival probability 
exponent for N — 2 is a = 1/2 independent of \x. We verify this result in our simulation. See Figure 5. Note that 
this result is very different from the result obtained in Ref. [21| where a depends on fi for N — 2 and higher. 
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We also simulate N = 3 vicious Levy flights. Because of the linearity of fractional derivatives, the mapping of 
two vicious Levy flights to a single Levy flights in an absorbing plane holds. However, due to the lack of rotational 
invariance of the Riesz-Feller derivative, the wedge mapping that holds for vicious walkers and now for vicious 
accelerating walkers, does not apply to vicious Levy flights. In order to make progress, since for N = 2 there exists 
a power-law distribution, we conjecture that the survival probability distribution scales as a power-law at long times 
for N > 2 and measure a. Figure 5 plots the survival probability exponents for N = 3 vicious Levy flights for several 
different Levy indices. For N — 2 all values of fj, yield the same survival probability exponent of 1/2, in agreement 
with the Sparre- Andersen theorem. However, the N = 3 exponents appear to vary with \i. For instance, for fi = 1, 
a = 1.31 ± 0.03. While the 0.19 difference between /i = 1.0 and \i = 2 is small, the difference grows with N. For 
example, for TV = 4 and fj, = 1, a = 2.3 ± 0.1 and for fj, = 2, a = 2.91 ± 0.09. Based on this data, we speculate that 
for N > 2, a depends on \i. 

A few comments on the technical aspects of the simulations are in order. We implement an upper cut-off on the Levy 
steps so that at long enough times, the survival probability distribution approaches the random walker result [24j. 




The convergence also depends on the Levy index. For example, for /i = 1 and stepsize cut-off S c = 10 — 100, the 
convergence to the random walker result is fast such that the asymptote to a power-law beyond t w 10 2 is in agreement 
with the random walker result to within one standard deviation. In contrast for ^ = 1.6, S c — 10 9 , and time scales 
beyond t sa 10 8 , convergence to the random walker result is observed. Secondly, for fi = 1 we also generated Cauchy 
distributed numbers directly and found good agreement with the power-law generated /j, = 1 result. 



To summarize, we have generalized the vicious walker problem in two different ways: (1) vicious accelerating walkers 
and (2) vicious Levy flights as defined by non-crossing. For both generalizations, the typical analytical technique of the 
method of images fails. Analytical results for N = 2 are readily obtainable since both problems can be mapped to the 
first passage problem of a single accelerating walker or Levy flight with the appropriate absorbing boundary or region. 
We demonstrate that the N vicious walker mapping to one walker in N — 1 dimensions in a wedge geometry generalizes 
to vicious accelerating walkers. We also conjecture, based on our numerical data, that there exists an upper bound 
on the survival probability exponent of a = |A(A — 1) for N vicious accelerating walkers. An analytical calculation 
for the N — 3 case corresponding to one acclerating walker in a two-dimensional wedge geometry would be the next 
logical step. The heuristic argument for the absorbing accelerating walker in one-dimension using a new time counter 
and Levy flights may eventually become useful to analyze the two-dimensional wedge problem, however, it is not 
yet clear how to generalize the Riesz-Feller derivative to dimensions higher than unity. There also exists a recent 
numerical result for the survival probability distribution for the two-dimensional Fractional Brownian motion process, 
originally introduced by Kolmogorov [25| , in a wedge [26} . We anticipate more study of non-Markovian processes in 
dimensions higher than unity in the near future. Indeed, a non-Markovian extension of Dyson's Brownian motion 
model to, for example, include inertia, may be related to N vicious accelerating walkers to arrive at a new class of 
random matrices. It may also be interesting to investigate other ordering problems of randomly accelerating walkers 
on a line such as the Gaussian equivalent of the "leader" and the "laggard" problem [27| . 

Finally, given our numerical results, we speculate that the survival probability exponent for N vicious Levy flights 
(as defined by no-crossing) depends on /i for fixed N > 2. We also refer to a new result where vicious Levy flights 
are defined as annihilating when any two Levy flights come within some range of each other and a depends on /i 
even for N = 2 (2lj . While the survival probability exponent in one-dimension is independent of the Levy index, 
as a consequence of the powerful Sparre- Andersen theorem, we anticipate that this superuniversality may be broken 
in dimensions higher than unity and the universality of each Levy index becomes exposed. In light of our results, a 
higher dimensional generalization (or modification) of the Sparre- Andersen theorem should be on the forefront of at 
least several statistical physicists and mathematicians minds. 

JMS would like to acknowledge M. Jeng for performing some preliminary simulations on no-crossing vicious Levy 
flights, Eli Hawkins for helpful discussion, and The Aspen Center of Physics where some of this work was conducted. 
JMS is supported by NSF-DMR-0645373. 
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